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Abstract 


A  methodology  has  been  developed  at  DRDC  Valcartier  for  the  calibration  of  infrared 
(IR)  cameras  when  taking  IR  signature  measurements  of  military  targets.  This 
methodology  is  the  basis  of  WinISAS  (Windows  Infrared  Signature  Analysis 
Software)  developed  at  DRDC  Valcartier  for  the  radiometric  analysis  of  IR  imagery 
and  consists  of  a  radiometric  calibration  of  the  grey  levels  in  the  image,  using  extended 
reference  blackbody  sources.  The  methodology  has  been  developed  and  used  mainly 
for  the  IR  signature  analysis  of  extended  targets,  fully  resolved  by  the  camera  IR 
sensor.  The  purpose  of  this  memorandum  is  to  demonstrate  that  the  same  methodology 
and  WinISAS  software  can  be  applied  to  the  measurement  of  point  targets,  which  are 
targets  with  an  angular  coverage  smaller  than  the  individual  detector  element  angular 
coverage  in  an  image  acquired  from  a  two-dimensional  IR  focal  plane  array  camera. 


Resume 


RDDC  Valcartier  a  mis  au  point  une  methodologie  pour  Tetalonnage  de  cameras  a 
infrarouge  (IR)  lors  de  la  prise  de  mesures  de  la  signature  IR  de  cibles  militaires.  Cette 
methodologie  est  la  base  du  logiciel  WinISAS  (Windows  Infrared  Signature  Analysis 
Software)  developpe  a  RDDC  Valcartier  pour  Tanalyse  radiometrique  d’images  IR  et 
consiste  en  un  etalonnage  radiometrique  des  niveaux  de  gris  de  T  image  a  partir  de 
corps  noirs  de  reference  de  grandes  dimensions.  Cette  methodologie  a  ete  mise  au 
point  principalement  pour  Tanalyse  de  la  signature  IR  de  cibles  etendues, 
completement  resolues  par  le  capteur  IR  de  la  camera.  L’objectif  de  ce  memorandum 
est  de  demontrer  que  la  meme  methodologie  ainsi  que  le  logiciel  WinISAS  peuvent 
s’appliquer  a  la  mesure  de  cibles  ponctuelles,  c’est-a-dire  des  cibles  de  couverture 
angulaire  plus  petite  que  celle  d’un  element  de  la  matrice,  dans  une  image  obtenue  par 
une  matrice  bi-dimensionnelle  de  detecteurs  IR  au  plan  focal  de  la  camera. 


DRDC  Valcartier  TM  2004-076 


I 


This  page  intentionally  left  blank. 


DRDC  Valcartier  TM  2004-076 


Executive  summary 


A  methodology  has  been  developed  at  DRDC  Valcartier  for  the  calibration  of  infrared 
(IR)  cameras  when  taking  IR  signature  measurements  of  military  targets.  This 
methodology,  documented  in  Refs.  1  and  2,  is  the  basis  of  WinISAS  (Windows 
Infrared  Signature  Analysis  Software),  developed  at  DRDC  Valcartier  for  the 
radiometric  analysis  of  IR  imagery.  It  consists  of  a  radiometric  calibration  of  the  grey 
levels  in  the  image,  using  extended  reference  blackbody  sources.  The  methodology  has 
been  developed  and  used  mainly  for  the  IR  signature  analysis  of  extended  targets,  fully 
resolved  by  the  camera  IR  sensor. 

The  Naval  Engineering  Test  Establishment  (NETE)  has  used  WinISAS  for  the  analysis 
of  the  IR  data  from  many  trials,  and  more  recently,  for  the  analysis  of  the  co-recorded 
calibrated  IR  video  data  collected  during  the  SIRIUS  trials  in  the  Bahamas  onboard  the 
HMCS  Montreal,  from  6  to  24  October  2003.  Questions  have  been  raised  about  the 
applicability  of  the  DRDC  Valcartier  calibration  and  analysis  methodology  when  using 
extended  (well-resolved)  reference  sources  to  calibrate  the  IR  cameras  and  obtain  the 
radiant  intensity  contrast  (A4)  of  a  point  target. 

The  purpose  of  this  memorandum  is  to  demonstrate  that  the  same  methodology  and 
WinISAS  can  be  applied  for  the  IR  analysis  of  extended  as  well  as  point  targets,  which 
are  targets  with  an  angular  coverage  smaller  than  the  individual  detector  element 
angular  coverage  in  an  image  acquired  from  a  two-dimensional  IR  focal  plane  array 
camera. 

Since  no  real  target  is  ever,  in  practice,  a  point  target  (because  it  would  need  to  be 
infinitesimally  small),  the  notion  of  a  point  target  is  always  defined  relatively  to  the 
detector  element’s  projected  area  at  the  target  range  and  the  projected  area  of  the 
imager’s  point-spread  function  (PSF).  The  intensity  of  any  target  is  defined  as  its 
radiance  times  its  cross-section  area.  By  using  an  extended-source  radiance  calibration 
method  to  evaluate  the  intensity  of  a  point  target,  it  has  been  shown  experimentally 
that,  for  a  point  target,  an  apparent  radiance  value  is  effectively  determined  instead  of 
the  real  target  radiance  because  the  latter  does  not  completely  fill  the  detector  element. 
To  this  apparent  radiance  is  associated  a  projected  detector  element  area  that  is  larger 
than  the  actual  target  area,  and  the  exact  radiant  intensity  value  is  obtained  when  the 
two  values  are  multiplied.  The  laboratory  measurements  also  demonstrate  that  the 
position  of  the  target  within  the  field  of  view  or  the  focus  setting  of  the  camera  have 
little  effect  on  the  calculated  intensity  values,  as  expected,  because  the  broadening  of 
the  target  image  caused  by  distortions  or  focus  settings  distributes  the  energy  over 
more  detector  elements,  thus  compensating  a  decrease  in  signal  by  a  larger  apparent 
target  area. 
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iii 


DRDC  Valcartier  TM  2004-076 


Sommaire 


RDDC  Valcartier  a  mis  au  point  une  methodologie  pour  I’etalonnage  de  cameras  a 
infrarouge  (IR)  lors  de  la  prise  de  mesures  de  la  signature  IR  de  cibles  militaires.  Cette 
methodologie  est  documentee  aux  ref  1  et  2  et  est  a  la  base  du  logiciel  WinISAS 
(Windows  Infrared  Signature  Analysis  Software),  egalement  mis  au  point  a  RDDC 
Valcartier  pour  I’analyse  radiometrique  d’images  IR.  Elle  permet  I’etablissement  d’une 
correspondance  radiometrique  avec  les  niveaux  de  gris  dans  I’image,  a  partir  de  corps 
noirs  de  reference  de  grandes  dimensions.  Cette  methodologie  a  ete  mise  au  point 
principalement  pour  1’ analyse  de  la  signature  IR  de  cibles  etendues,  completement 
resolues  par  le  capteur  IR  de  la  camera. 

Le  Centre  d’essais  techniques  (Mer)  (CETM)  a  utilise  WinISAS  pour  I’analyse  des 
donnees  IR  de  plusieurs  essais,  et  plus  recemment  pour  I’analyse  des  donnees  video  IR 
co-enregistrees  lors  des  essais  du  systeme  SIRIUS  dans  les  Bahamas,  a  bord  du  NCSM 
Montreal,  du  6  au  24  octobre  2003.  La  methodologie  d’etalonnage  et  d’analyse  de 
RDDC  Valcartier,  basee  sur  I’utilisation  de  sources  de  reference  etendues  (resolues) 
pour  I’etalonnage  des  cameras  IR,  a  ete  mise  en  question  quant  a  son  applicabilite  au 
calcul  du  contraste  en  intensite  (Alt)  de  cibles  ponctuelles. 

Ce  memorandum  a  pour  objectif  de  demontrer  que  la  meme  methodologie  ainsi  que  le 
logiciel  WinISAS  peuvent  etre  appliques  tant  a  I’analyse  de  cibles  etendues  qu’a  des 
cibles  ponctuelles,  c’est-a-dire  des  cibles  dont  les  dimensions  angulaires  sont  plus 
petites  que  cedes  de  1’ element  individuel  d’une  matrice  de  detecteurs  IR  au  plan  focal 
de  la  camera. 

Puisqu’en  pratique,  aucune  cible  n’est  vraiment  ponctuelle  (il  faudrait  qu’elle  soit 
infiniment  petite),  la  notion  de  cible  ponctuelle  est  toujours  definie  relativement  a  la 
surface  projetee  d’un  element  de  la  matrice  de  detecteurs  a  la  distance  de  la  cible  ou  a 
la  surface  projetee  de  la fonction  d’etalement  de  la  camera.  L’intensite  de  n’importe 
quelle  cible  est  definie  comme  sa  radiance  energetique  multipliee  par  sa  surface.  On 
demontre  experimentalement  qu’en  utilisant  une  methode  d’etalonnage  en  radiance  au 
moyen  de  sources  etendues  dans  le  cas  d’une  cible  ponctuelle,  on  obtient  une  valeur  de 
radiance  apparente  au  lieu  de  la  radiance  reelle  de  la  cible  parce  que  I’image  de  celle- 
ci  ne  remplit  pas  un  element  de  la  matrice  de  detecteurs.  A  cette  radiance  apparente, 
on  associe  la  surface  projetee  d’un  element  du  detecteur  qui  est  plus  grande  que  la 
dimension  reelle  de  la  cible  et  on  obtient  la  valeur  exacte  de  1’ intensite  radiometrique 
lorsque  les  deux  valeurs  sont  multipliees.  Des  mesures  en  laboratoire  demontrent  aussi 
que  la  position  de  la  cible  dans  le  champ  de  vision  de  la  camera  ainsi  que  I’ajustement 
de  mise  au  foyer  ont  peu  d’effet  sur  la  valeur  calculee  de  I’intensite,  puisque 
I’elargissement  de  I’image  de  la  cible  cause  par  des  distorsions  ou  une  mauvaise  mise 
au  point  distribue  I’energie  sur  plusieurs  elements  detecteurs,  compensant  de  ce  fait  la 
diminution  du  signal  par  une  plus  grande  surface  apparente  de  la  cible. 
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1.  Introduction 


A  methodology  has  been  developed  at  DRDC  Valcartier  for  the  calibration  of  infrared 
(IR)  cameras  when  taking  IR  signature  measurements  of  military  targets.  This 
methodology  is  documented  in  Refs.  1  and  2  and  is  the  basis  of  WinISAS  (Windows 
Infrared  Signature  Analysis  Software),  developed  at  DRDC  Valcartier  for  the 
radiometric  analysis  of  IR  imagery.  It  consists  of  a  radiometric  calibration  of  the  grey 
levels  in  the  image,  using  extended  reference  blackbody  sources.  The  methodology  has 
been  developed  and  used  mainly  for  the  IR  signature  analysis  of  extended  targets,  fully 
resolved  by  the  camera  IR  sensor. 

The  Naval  Engineering  Test  Establishment  (NETE)  has  used  WinISAS  for  the  analysis 
of  the  results  of  many  trials,  and  more  recently,  for  the  analysis  of  the  co-recorded 
calibrated  video  data  collected  during  the  SIRIUS  trials  in  the  Bahamas  onboard  the 
HMCS  Montreal,  from  6  to  24  October  2003.  During  these  trials,  suspended  targets 
were  towed  by  a  helicopter,  and  a  fighter  aircraft  was  also  used  to  test  the  SIRIUS 
system.  The  suspended  targets  were  considered  to  be  point  targets  at  longer  ranges, 
and  the  target  fighter  aircraft  can  also  be  considered  to  be  a  point  target  at  its 
maximum  detection  range  and  for  a  good  part  of  its  trajectory.  Questions  have  been 
raised  about  the  applicability  of  the  DRDC  calibration  and  analysis  methodology  when 
using  extended  (well  resolved)  reference  sources  to  calibrate  the  cameras  and  obtain 
the  radiant  intensity  contrast  (Alt)  of  a  point  target. 

The  purpose  of  this  memorandum  is  to  demonstrate  that  the  same  methodology  and 
WinISAS  can  be  applied  to  the  measurement  of  extended  as  well  as  point  targets, 
which  are  targets  with  an  angular  coverage  smaller  than  the  individual  detector 
element  angular  coverage  in  an  image  acquired  from  a  two-dimensional  IR  focal  plane 
array  camera. 

2.  Basic  relations 


The  reader  is  referred  to  Refs.  1  and  2  for  a  complete  description  of  DRDC  Valcartier 
extended  source  calibration  methodology.  In  brief,  two  external  blackbody  reference 
sources  of  known  temperatures  are  placed  at  some  known  distance  from  the  cameras 
and  images  of  these  sources  are  recorded  prior  to  or  immediately  after  recording 
images  of  a  target.  The  sources  must  be  well  resolved  on  the  detector  so  that  its  solid 
angle,  as  seen  by  the  cameras,  covers  several  detector  element  solid  angles,  ^^/(see 
Fig.  1).  Whenever  possible,  if  the  reference  sources  are  larger  than  the  camera’s  lens 
diameter,  they  can  also  be  placed  directly  against  the  lens  for  calibration  purpose 
(calibration  source  covers  all  the  detector  elements).  Furthermore,  the  gain  and  offset 
settings  of  the  camera  used  during  the  calibration  should  also  be  used  when  recording 
the  target,  if  there  is  no  way  of  converting  the  calibration  from  one  setting  to  another. 
The  basic  idea  in  calibration  is  to  establish  a  one-to-one  correspondence  between  the 
apparent  radiance  L  (W/m^-sr)  from  any  part  of  the  imaged  scene  and  the  grey  (signal) 
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level  measured  by  the  camera,  assuming  that  the  detector  responds  linearly  to  the 
incident  radiometric  flux  (W),  which  in  turn  is  proportional  to  the  target  radiance 
integrated  over  the  limits  of  the  camera’s  spectral  response  R(A).  A  target  apparent 
radiance  is  the  radiance  measured  by  the  camera  for  a  target  at  a  known  range  and 
including  the  effect  of  atmospheric  attenuation.  WinISAS  is  designed  to  take  the 
amount  of  atmospheric  attenuation  into  account  and  to  correct  for  its  effect  to  obtain 
the  unattenuated  radiance  of  the  target.  Of  course,  when  an  IR  image  analysis  is 
performed  to  extract  the  IR  signature  of  a  target,  in  most  cases  the  results  are  given  in 
contrast  values,  i.e.  the  apparent  or  corrected  radiance  of  the  background  (selected  in 
the  scene  immediately  surrounding  the  target)  is  evaluated  and  subtracted  from  the 
measured  apparent  or  corrected  radiance  of  each  target  pixel.  At  this  point,  for 
simplicity,  the  atmospheric  attenuation  effect  is  ignored  when  describing  the  basic 
relations  between  the  target  radiance  Lt,  the  incident  flux  Oto  and  the  irradiance  Eto  on 
the  camera  optics  from  the  target,  and  the  flux  Ot^  on  a  detector  element,  when  the 
extent  of  the  target  is  much  larger  than  Q^.  The  exact  relationships  are  given  by  Eqs 
2.1  and  2.2,  and  their  geometry  is  illustrated  in  Fig.  1. 


o  = 

^  to 

f  O,  ^ 

E,  = 

to 

.  A 

\  0  J 

where: 

\R  j 


r  A  \ 


=  L. 


\R^j 


(2.1) 

(2.2) 


Lt  is  the  target  radiance  integrated  over  the  spectral  response  of  the  camera; 

Adt  is  the  area  of  the  sensitive  surface  of  the  detector  element  projected  on  the 
target; 

Ao  is  the  area  of  the  camera  front  optics; 

R  is  the  range  of  the  target  from  the  camera 

To  obtain  contrast  values,  we  subtract  the  background  corresponding  signal  to  the 
target  signal: 

A(D,,  =(L,  (2.1b) 

AE,^=AL,Q,  (2.2b) 
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Figure  1.  Geometry  describing  the  relationship  between  target  and  camera  parameters 

Assuming  no  loss  and  no  degradation  of  the  image  by  the  optics  (geometric  image), 
the  radiometric  flux  contrast  AOt^  on  a  detector  element  from  the  extended  target  is  the 
same  as  AOto  on  the  front  optics,  originating  from  the  detector  element’s  projected  area 
in  the  target  plane.  It  is  obtained  from  the  following  relation: 


A<D,^  =  = 


rAi,4Z)0 

n 

f  4 

I  f  ) 

\  J 

(2.3) 


where: 

Ad  is  the  detector’s  sensitive  area  in  the  focal  plane  of  the  optics; 

/  is  the  focal  length  of  the  optics; 

D  is  the  diameter  of  the  optics; 

F#  is  the  optics  f-number  if/D)  of  the  optics. 

From  Eq.  2.3  it  is  obvious  that  is  directly  proportional  to  the  target  radiance  Lt,  as 
the  other  parameters  are  camera  speciflc  constants. 

3.  Extended  target  and  effect  of  the  PSF 

For  an  exact  analysis  of  the  radiometric  flux  hitting  a  detector  element,  the  effect  of 
image  degradation  caused  by  the  optical  point-spread  function  {PSF)  must  be  taken 
into  account.  In  any  imaging  system,  the  observed  image  results  from  the  convolution 
of  the  perfect  geometric  image  of  the  scene  with  the  PSF  of  the  front  optics,  also 
known  as  the  blur  spot.  This  is  the  best  focal  point  that  can  be  obtained  when  the 
effects  of  diffraction  and  optical  aberration  are  considered.  For  simplicity,  during  the 
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following  discussions,  we  will  assume  that  the  transmission  factor  of  the  front  optics  is 

100%. 


The  following  equation  expresses  the  flux  contrast  received  by  a  detector  element  from 
an  extended  target  when  the  effect  of  the  PSF  is  introduced  into  equation  2.3: 


TT 

4T 


][PSFix,y)*M,ix,y)\  ec^l 


^  X  ^ 


rect 


(y\ 

yLj 


dxdy 


(3.1) 


where: 

ALti(x,y)  is  the  image  mapping  of  the  target  radiance  contrast  in  the  focal  plane; 

Xd,  Yd  are  the  respective  horizontal  and  vertical  dimensions  of  the  detector 
element’s  sensitive  area; 


rect 


-  =‘ 

X) 


for 


- <x<  — 

2  2  0  otherwise 


and  the  convolution  (*  operator)  of  the  target  radiance  image  with  the  PSF  is  given  by: 
PSF{x,  y)  *  (x,  ;^)  =  j  PSF{^  -  x,  77  -  y )  x  (^,  77)6? rj  (3.2) 


L 


Figure  2.  Representation  of  steps  in  equation  3.1 
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Equation  3.2  shows  that  the  convolution  is  performed  over  the  entire  target  image 
Lti(x,y)  (calibrated  in  target  radiance  value).  However,  when  the  target  is  larger  than  the 
detector  area,  then  the  integration  of  the  convolution  function  in  Eq.  3.1  is  limited  to 
the  detector’s  dimension  by  the  rect( )  functions.  Thus,  Eq.  3.1  represents  the  total  flux 
received  by  a  detector  element  from  a  uniform  target  that  is  larger  in  size  than  the  solid 
angle  seen  by  the  detector  element.  This  value  corresponds  to  the  area  under  the 
portion  of  the  convolution  curve  delimited  by  the  detector’s  sensitive  area,  as  shown  in 
Figure  2. 

4.  Point  target  and  the  effect  of  the  PSF 

. sJf . 


The  same  equation,  Eq.  3.1,  applies  to  the  description  of  the  total  flux  contrast 
resulting  from  a  point  target  on  a  detector  element.  However,  in  this  case  the  target 
image  radiance  contrast  function  lsLti(x,y)  is  signiflcantly  smaller  than  the  detector 
size,  as  shown  in  Figure  3.  Then  the  convolution  of  the  target  image  lSLti(x,y)  with  the 
PSF(x,y)  yields  a  function  that  is  very  similar  to  the  PSF(x,y)  function,  except  for  a 
slightly  larger  width. 


Thus,  the  total  flux  contrast  on  a  detector  element  is  the  integrated  volume  under 
the  convolution  surface  (when  considering  the  two  dimensions).  As  the  spatial  extent 
of  the  convolution  is  smaller  than  the  detector’s  area  (when  the  image  of  the  target  is 
centered  on  a  detector  element),  the  measured  flux  is  not  affected  by  the  size  of  the 
detector  element,  as  opposed  to  the  case  of  an  extended  target.  Consequently,  when  the 
detector  has  been  calibrated  using  an  extended  reference  source  of  constant  radiance, 
the  measured  flux  contrast  from  a  point  target  that  falls  entirely  within  one  detector 
element  can  be  interpreted  as  that  coming  from  an  extended  target  of  constant  apparent 
radiance  contrast  (ALt_app)  that  is  proportional  to  the  measured  flux  contrast,  AO^^, 
divided  by  the  detector’s  sensitive  area,  similarly  to  Eq.  2.3.  In  this  case,  the  apparent 
radiance  contrast  is  given  by 


AL. 


't_app 


n 


2  Y 


AO. 


\  j 


(4.1) 


and  will  always  be  less  than  the  true  radiance  contrast,  AL^,  of  the  target.We  will  see  in 
the  next  chapter  that  this  lower  apparent  radiance  value,  multiplied  by  a  larger 
projected  area  than  the  actual  target  size,  can  give  the  real  target  radiant  intensity 
contrast  value  A/^  (W/sr). 


A/,  =  AL, 

t  t  _  opp  ut 


(4.2) 
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X 


X 


Equivalent  to  averaged  over  detector  area 


Figure  3.  Representation  of  steps  in  equation  3. 1  for  a  point  target 


5.  Calibration  with  an  extended  reference  source  and 
intensity  measurements 


5.1  Case  of  large  target 

Assuming  that  the  measured  signal  from  each  detector  element  (grey  level,  g)  is 
proportional  to  the  flux  contrast,  AO^^,  incident  on  the  detector’s  sensitive  area,  and 
that  the  calibration  has  been  performed  using  an  extended  source,  then  the  measured 
flux  can  be  directly  related  to  the  known  radiance  of  the  reference  source  within  the 
spectral  response  of  the  camera  (Refs.  1  and  2). 
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At  this  point,  if  the  camera  has  a  2D  focal  plane  detector  array,  the  notion  of fill  factor 
must  be  introduced.  In  most  detector  arrays,  the  detector’s  sensitive  area  does  not 
cover  all  the  space  between  adjacent  detector  elements  because  a  part  of  that  space  is 
used  for  the  layout  of  electrodes  or  wires  to  connect  the  photosensitive  part  to  the 
electronic  readout  circuits,  thus  creating  some  dead  space  around  the  sensitive  area  of 
each  detector  element.  The  ratio  between  the  actual  sensitive  area  of  the  detector  array 
and  the  total  area  of  the  detector  array  is  called  the  fill  factor  (also  the  ratio  of  the 
sensitive  area  of  a  detector  element  to  the  total  space  available  for  each  element). 

When  the  camera  outputs  a  direct  digital  image  from  a  2D  focal-plane  array,  the 
following  procedure  can  be  used  to  obtain  the  overall  intensity  contrast,  A/^,  of  a 
uniform  target  whose  extent  is  much  larger  than  the  solid  angle  coved  by  a  single 
detector  element. 

Using  the  radiance  calibration  of  the  camera  by  means  of  an  extended  source,  the 
overall  target  intensity  is  calculated  from  the  summation  of  the  corresponding  radiance 
values  for  all  the  detector  elements  on  target,  multiplied  by  their  projected  area  in  the 
target  plane  with  a  fill  factor  of  one  (Adtj).  This  is  better  accomplished  by  summing  the 
histogram  of  the  detector  elements  on  target  as  a  function  of  their  grey  levels,  Hfg), 
multiplied  by  the  corresponding  radiance,  L(g)  and  detector  projected  area^^^/(Eq. 
5.1).  The  use  of  a  fill  factor  of  one  when  projecting  the  detector  area  onto  the  target 
plane  is  based  on  the  premise  that  the  target  radiance  value  is  the  same  in  the  dead 
space  between  the  detector  elements  as  it  is  on  the  detector  element  sensitive  area.  The 
overall  target  radiant  intensity  f  (W/sr)  is  then  obtained  from 

I,=Y,L{g)HXg)A,,^  (5.1) 


where: 

L(g)  is  the  calibrated  radiance  corresponding  to  grey  level  g; 

Hfg)  is  the  histogram  of  the  detector  elements  on  target  as  a  function  of  grey  level 
(number  of  detector  elements  with  grey  level  g); 

the  projected  area  of  the  detector  element  in  the  target  plane,  with  a  fill  factor 

of  1. 

The  target  radiant  intensity  contrast  is  obtained  by  subtracting  an  equivalent  area  at  the 
background  radiance  value: 


A/,  =  -  'LhH,ig)Af  = 

g  g  g 

AI,=X^(g)H,(g)A,^^ 

g 


(5.2) 
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5.2  Case  of  point  target  and  PSF  smaller  than  detector 
sensitive  area 


A  target  is  considered  as  a  point  target  when  its  extension  is  smaller  than  the  projected 
PSF  of  the  optical  system  and,  in  most  cases  also  smaller  than  the  detector  element’s 
pro jQctQd  sensitive  area.  In  that  case,  Eq.  5.1  becomes: 

(5-3) 

where: 

Lt_app  is  the  apparent  radiance  measured  by  the  detector  element,  based  on  its 
calibration  with  an  extended  source; 

Adt  is  the  projected  area  of  the  sensitive  part  of  the  detector  element  in  the  target 
plane. 

The  contrast  radiant  intensity  is  obtained  similarly  to  Eq.  5.2: 

-  m_app  ~  ^b)^dt  -  ^t_app^dt  (5-4) 

Since,  upon  calibration  with  an  extended  reference  source,  the  flux  was  integrated  over 
the  detector’s  sensitive  area,  and  since  the  point  target  is  smaller  than  that  area,  it  is 
normal  that,  to  recover  the  radiant  intensity  It  of  the  target,  the  projected  sensitive  area 
of  the  detector  element  is  used.  Thus,  the  apparent  radiance  value  Lt_app  is  lower  than 
the  target’s  real  radiance  value  and  is  compensated  by  multiplying  with  a  projected 
area  larger  than  the  actual  target  area  to  obtain  an  exact  target  radiant  intensity  value. 

Note:  This  case  is  almost  hypothetical  because  an  optical  designer  will  normally 
correct  his  design  until  the  PSF  approximately  matches  the  dimensions  of  the  detector 
element.  Bettering  the  optics  beyond  this  limit  only  increases  costs  without  any 
additional  benefits.  Moreover,  beyond  some  relatively  short  ranges,  the  turbulence 
effects  of  the  atmosphere  will  broaden  the  point  target  image  to  dimensions  exceeding 
the  span  of  a  single  detector  element. 

5.3  Case  of  point  target  and  PSF  larger  than  detector 
sensitive  area 

In  the  case  when  the  PSF  of  the  optics,  or  the  PSF  of  the  optics  combined  with  the 
atmospheric  turbulence  dispersion  effect,  is  larger  than  one  detector  element,  then  Eq. 
5.1  also  applies  to  point  targets  since  the  PSF  then  overflows  on  neighboring  detector 
elements.  Consequently,  the  summation  is  performed  over  all  the  detector  elements 
included  in  the  area  covered  by  the  PSF  convolved  with  the  point  target  (atmospheric 
turbulence  effect  included),  and  is  used  (detector  area  with  a  fill  factor  of  1.0)  to 
compensate  for  the  loss  of  energy  falling  on  the  dead  space  between  detector  elements. 

However,  in  such  a  case,  an  effective  fill  factor  smaller  than  one  may  have  an  impact 
on  the  accuracy  of  the  measured  intensity  value  if  only  a  few  pixels  are  included  in  the 
PSF.  Figure  4  illustrates  this  effect. 
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collected  energy  on  detector  elements  with  real  fill  factor  <1.0 

pm _ 

X 


L 


estimated  intensity  (assuming  a  fill  faetor  =  1.0) 


Figure  4.  Effect  of  an  effective  fill  factor  smaller  than  one  on  the  estimated  intensity 

As  can  be  seen,  the  detector  elements  are  collecting  radiant  energy  only  over  their 
active  area,  if  the  fill  factor  is  smaller  than  one  and  the  detector  element  integrating  on 
the  edge  of  the  PSF,  where  there  is  a  strong  gradient,  the  measured  radiance 
corresponds  to  the  average  radiance  over  the  actual  detector  collecting  area.  When  this 
apparent  radiance  is  converted  to  a  radiant  intensity  value,  the  average  radiance  is 
multiplied  by  the  projected  detector  area  assuming  a  fill  factor  of  one,  therefore  the 
same  measured  apparent  radiance  is  given  to  the  gap  adjacent  to  the  detector  element. 
Depending  of  the  position  of  the  detector  element  with  respect  to  the  gradient  signal, 
this  leads  to  either  an  overestimate  or  an  underestimate  of  the  true  intensity  value,  the 
smaller  the  number  of  pixel  elements  encompassed  in  the  PSF,  the  larger  the  error.  The 
error  will  also  vary  with  shifting  positions  of  the  target  over  the  detector  matrix.  As  the 
number  of  included  elements  increases  (more  samples),  this  edge  effect  tends  to 
becomes  negligible. 


5.4  Case  of  digitized  analog  image 

Oftentimes,  the  digital  signal  is  converted  to  an  NTSC  (or  PAL)  analog  video  output 
that  is  recorded  onto  analog  tape  before  being  converted  back  to  digital  format  using  a 
frame-grabber  card  for  further  analysis.  In  this  case,  no  spatial  relationship  exists 
between  the  digital  pixels  and  the  actual  detector  elements  in  the  focal-plane  array. 
Consequently,  calibration  with  an  extended  reference  source  should  be  performed 
following  the  digitization  of  the  images  with  the  frame-grabber  such  that  the  electronic 
bandpass  of  the  analog  video  amplifier,  the  recorder  and  the  frame-grabber  are 
included.  The  above  equations  may  then  be  applied  using  the  projected  area  of  the 
pixels  defined  by  the  frame-grabber  instead  of  the  projected  area  of  the  detector 
elements.  Most  likely,  the  additional  electronic  bandpass  will  have  smoothed  the 
discrete  signal  of  the  detector  elements  to  a  continuous  signal  during  the  digital  to 
analog  conversion  so  that,  in  this  case,  the  fill  factor  is  100%  and  makes  the  individual 
elements  indistinguishable.  However,  caution  should  always  be  exercised  to  make  sure 
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that  this  is  actually  the  case.  If  not,  for  point  targets,  an  effective  fill  factor  should  be 
determined  from  the  analog  and  digital  signals. 

6.  Experimental results 

To  verify  the  above  described  principles,  and  the  applicability  of  the  calibration 
methodology  and  WinISAS  to  point  target  measurements,  experimental  measurements 
were  taken  in  the  laboratory  with  known  blackbody  references  and  point  targets. 
Figure  5  illustrates  the  experimental  set-up. 


reference  2 


reference  1 


point  target 
pin  hole 


Figure  5.  Experimental  setup  for  calibration  and  point  target  measurements 

The  experimental  measurements  were  performed  with  an  LWIR  Phoenix  camera  from 
Indigo,  a  focal  plane  array  of  quantum  well  detectors,  640  x  512  elements,  operating  in 
the  8.0  -  9.2  waveband,  with  a  14-bit  digital  output,  using  an  100-mm  objective  lens 
which  gave  a  horizontal  field  of  view  of  9.2  degrees.  Two  extended  reference 
blackbodies  at  known  temperatures  were  used  for  the  radiometric  calibration  of  the 
camera.  The  target  was  a  velvet  black  painted  plate  in  which  discs  with  pin  holes  of 
various  size  could  be  inserted,  in  front  of  a  controlled  blackbody  whose  temperature 
could  be  adjusted  to  obtain  a  reasonable  signal  to  noise  ratio.  The  dimensions  of  the 
plate  holding  the  target  discs  could  also  be  used  for  the  spatial  calibration,  that  is  to 
obtain  the  pixel  angular  coverage  and  projected  area,  knowing  the  distance  of  the  plate 
(3.0  meters)  to  the  camera.  One  of  the  reference  blackbodies  was  kept  at  ambient 
temperature  (25°C),  while  the  other  was  set  at  90°C  for  an  optimum  coverage  of  the 
camera  dynamic  range.  Following  the  grabbing  of  calibration  images  (reference 
blackbodies  and  plate  of  known  dimensions),  target  images  were  grabbed  for  pinholes 
of  different  sizes,  from  fully  resolved  to  sub-pixel  coverage.  For  the  sub-pixel  targets, 
images  were  grabbed  with  the  target  appearing  at  different  positions  in  the  field  of 
view,  and  with  different  focus  settings,  in  order  to  see  the  effect  of  lens  distortions  and 
focus  adjustment. 
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Figure  6  Point  target,  0.17  mrad  diameter,  centered  and  in-focus 


Figure  7  Point  target,  0.17  mrad  diameter,  centered  and  slightly  out  of  focus 
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Figure  8  Point  target,  0.17  mrad  diameter,  to  the  right  and  with  coma  distortion 


Figure  9  Point  target,  0.17  mrad  diameter,  top  and  slightly  out  of  focus 
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Table  1  Real  and  measured  and  contrast  intensities  for  resolved  and  unresolved  targets,  using  WinISAS 
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Tgt  diam 

Tgt  IFOV 
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TgtL 

Bkgnd  L 

Tgt  area  (real)  Tgt  area  (meas]  Tgt  AL  (meas) 

Ag 
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A 
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6 
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90 
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NATO 

center 

6 

5.0800 
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90 

26.96 
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2.027E-05 

1.754E-04 

2.1293 
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3.735E-04 

3.668E-04 

1.81% 

NATO 

center 

6 

2.5400 

0.8621 

90 

26.96 

8.9382 

5.067E-06 

1.255E-04 

0.7519 
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9.436E-05 

9.132E-05 

3.34% 

NATO 

center 

6 

1.2700 

0.4310 

90 

26.96 

9.0149 

1.267E-06 

6.091  E-05 

0.3784 
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2.305E-05 

2.273E-05 

1 .39% 

NATO 

center,  weak  signal 

I  6 

0.7620 

0.2586 

200 

83.64 

8.9213 

4.560E-07 

7.223E-05 

0.4161 
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3.005E-05 

3.407E-05 

-11.79% 

NATO 

center,  weak  signal 

6 

0.7620 

0.2586 
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83.64 

8.9285 

4.560E-07 

7.223E-05 

0.4412 
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3.187E-05 

3.407E-05 

-6.46% 

NATO 

center,  weak  signal,  out  c 

I  6,1 

0.5080 

0.1724 
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293.27 
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2.027E-07 

5.565E-05 

0.885 
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4.925E-05 

5.760E-05 

-14.49% 

NATO 
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6.2 

0.5080 

0.1724 
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293.27 

9.0985 

2.027E-07 

6.623E-05 

0.8743 
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5.790E-05 

5.760E-05 

0.54% 

NATO 

center  right 

6,3 

0.5080 

0.1724 
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293.27 

9.0737 

2.027E-07 

6.071  E-05 

1.0197 
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6.1 91  E-05 

5.760E-05 

7.47% 
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right 

6,4 

0.5080 

0.1724 
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2.027E-07 

5.059E-05 
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5.232E-05 

5.760E-05 

-9.17% 

NATO 

center 

6,  5 

0.5080 

0.1724 

425 

293.27 

9.1225 

2.027E-07 

7.772E-05 

0.7625 
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5.926E-05 

5.759E-05 

2.90% 

NATO 

center  left 

6,6 

0.5080 

0.1724 

425 

293.27 

9.1110 

2.027E-07 

9.658E-05 

0.6513 

247 

6.290E-05 

5.759E-05 

9.22% 

NATO 

left 

6,  7 

0.5080 

0.1724 
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293.27 

9.1239 

2.027E-07 

6.071  E-05 

0.9581 
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5.817E-05 

5.759E-05 

1.00% 

NATO 

top 

6,  8 

0.5080 

0.1724 

425 

293.27 

9.0718 

2.027E-07 

7.772E-05 

0.7833 

363 

6.088E-05 

5.760E-05 

5.69% 

NATO 

bottom 

I  679 

0.5080 

0.1724 

425 

293.27 

9.1117 

2.027E-07 

7.772E-05 

0.6286 

396 

4.885E-05 

5.759E-05 

-15.17% 

NATO 

center,  on  pixel  ] 

6.  1 

0.5080 

0.1724 

425 

293.27 

9.1047 

2.027E-07 

7.175E-05 

0.7749 

293 

5.560E-05 

5.760E-05 

-3.47% 

NATO 

center,  out  of  focus 

6,2 

0.5080 

0.1724 

425 

293.27 

9.1012 

2.027E-07 

8.370E-05 

0.6632 

252 

5.551  E-05 

5.760E-05 

-3.62% 

NATO 

center  right,  out  of  focus 

6,  3 

0.5080 

0.1724 

425 

293.27 

9.0646 

2.027E-07 

6.477E-05 

0.8945 

339 

5.794E-05 

5.760E-05 

0.58% 

NATO 

right,  out  of  focus 

6,4 

0.5080 

0.1724 

425 

293.27 

9.1023 

2.027E-07 

7.772E-05 

0.7265 

276 

5.646E-05 

5.760E-05 

-1 .97% 

NATO 

center  left,  out  of  focus 

6,  5 

0.5080 

0.1724 

425 

293.27 

9.1189 

2.027E-07 

8.278E-05 

0.6867 

260 

5.685E-05 

5.759E-05 

-1 .30% 

NATO 

left,  out  of  focus 
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NATO 

center 

7,  1 

0.5080 

0.1724 

372 

234.58 

9.1767 

2.027E-07 

7.798E-05 

0.5354 

208 

4.175E-05 

4.568E-05 

-8.61% 

NATO 

center 

7,  2 

0.5080 

0.1724 

373 

235.64 

9.1710 

2.027E-07 

7.752E-05 

0.5046 

196 

3.912E-05 

4.590E-05 

-14.78% 
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Figures  6-9  illustrate  IR  images  of  the  point  target  at  different  locations  in  the  field  of 
view  and  with  slightly  different  focus  settings.  The  out-of-focus  cases  correspond  to  a 
signal  drop  by  a  factor  of  2.  In  Fig.  8,  a  slight  blurring  can  be  seen  caused  by  coma 
distortion. 

Table  1  lists  the  results  for  resolved  and  unresolved  targets  as  measured  in  the 
laboratory  with  the  setup  described  in  Fig.  5.  Column  3  gives  the  target  diameter  in 
mrad.  Each  detector  element  covers  an  instantaneous  field  of  view  of  0.25  mrad.  The 
temperature  of  the  blackbody  source  behind  the  target  plate  is  given  in  column  4.  The 
measured  and  calculated  (real)  target  intensity  contrasts  appear  in  columns  1 1  and  12, 
respectively,  while  the  percentage  of  error  is  shown  in  column  13.  The  yellow  shaded 
lines  indicate  errors  above  10%.  Except  for  a  case  where  the  signal  was  weak,  it  can  be 
observed  that  those  cases  occur  mainly  for  target  in- focus,  centered  on  a  pixel  and  in 
the  center  of  the  field  of  view,  which  can  be  explained  by  the  fill-factor  effect,  as 
shown  in  Fig.  4. 

Overall,  the  results  are  acceptable,  considering  all  the  potential  sources  of  error  such  as 
the  emissivity  of  the  calibration  sources  and  the  blackbody  source  behind  the  target 
plate  (assumed  to  be  1.0),  the  temperature  reading  of  the  blackbody,  the  atmospheric 
effect  in  the  optical  path  (assumed  negligible,  but  the  spectral  response  of  the  camera 
has  a  long  tail  down  to  6  pm),  the  digital  noise  of  the  camera,  the  fill  factor  effect, 
which  may  well  sum  up  to  ±  10%.  It  is  seen  that  the  position  of  the  target  within  the 
field  of  view  or  the  focus  setting  of  the  camera  have  little  effect  on  the  calculated 
intensity  values  as  expected,  because  the  broadening  of  the  target  image  caused  by 
distortions  or  focus  settings  distributes  the  energy  over  more  detector  elements,  thus 
compensating  for  a  decrease  in  the  apparent  measured  radiance  values  by  a  larger 
apparent  target  area. 

7.  Conclusion 


It  has  been  demonstrated  mathematically  that  one  may  use  an  infrared  camera  that  has 
been  radiometrically  calibrated  in  radiance  units  with  extended  reference  sources,  to 
measure  the  radiant  intensity  of  point  targets.  Since  no  real  target  is  ever,  in  practice,  a 
point  target  (because  it  would  need  to  be  infinitesimally  small),  the  notion  of  a  point 
target  is  always  defined  relative  to  the  detector  element’s  projected  area  at  the  target 
range  and  the  projected  area  of  the  imager’s  PSF.  The  intensity  of  any  target  is  defined 
as  its  radiance  times  its  cross-section  area.  By  using  an  extended-source  radiance 
calibration  method  to  evaluate  the  intensity  of  a  point  target,  it  has  been  shown  that  an 
apparent  radiance  value  is  effectively  determined  instead  of  the  real  target  radiance,  as 
the  latter  does  not  completely  fill  the  detector  element.  To  this  apparent  radiance  is 
associated  a  projected  detector  element  area  that  is  larger  than  the  actual  target  area, 
and  the  exact  radiant  intensity  value  is  obtained  when  the  two  values  are  multiplied. 
However,  some  caution  must  still  be  exercised  to  properly  perform  the  calibration 
procedure,  as  described  in  the  calibration  methodology  detailed  in  Refs.  1  and  2,  to 
avoid  saturation  of  the  reference  sources  signals,  as  well  as  of  the  target  signal  during 
the  recording.  When  a  camera  with  a  focal  plane  array  of  sensors  is  used  with  an 


14 


DRDC  Valcartier  TM  2004-076 


analog  recorder,  the  direct  relationship  between  the  true  detector  elements  and  the 
digital  pixels,  obtained  upon  digital  re-conversion  prior  to  analysis,  is  lost.  In  this  case, 
further  caution  should  also  be  exercised  to  ascertain  that  a  smooth  signal  is  obtained  in 
the  analog  video  when  scanning  a  point  source  across  the  focal  plane  array.  This  will 
prevent  the  effect  of  uncorrected  dead  spaces  between  detector  elements  when  using 
the  frame  grabber  re-converted  digital  pixels  for  analysis  in  lieu  of  the  true  detector 
digital  signal,  by  ensuring  that  the  PSF  is  actually  larger  than  a  single  detector’s 
sensitive  area  or  that  spurious  dead  space  effects  are  eliminated  by  the  electronics 
bandwidth. 

The  methodology  described  in  this  memorandum  has  been  used  successfully  by  NETE 
for  the  analysis  of  the  co-recording  measurements  taken  during  the  SIRIUS  trials  in 
October  2003  and  the  comparison  with  the  results  obtained  from  the  SIRIUS  system 
based  on  modeling  and  detector  performance  parameters  has  shown  a  more  than 
satisfactory  agreement^ 


^  D.  Dion,  Y.  Bedard,  L.  Forand,  E.  Fortier,  P.  Chevrette,  “DRDC  Analysis  of  the  Sirius  Warm 
Weather  Trial  Data”,  DRDC  Valcartier  report  to  be  published. 
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